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(54) Verfahren zur Verarbeitung eines Eingangsbiides 

(57) Die Erfindung bezieht sich auf ein Bildverarbei- 
tungsverlahren bei dem ein Eingangsbild zun^chst in 
eine Reihe von Detailbildern zerlegt wird. denen jeweils 
ein Ortsfrequenzbereich zugeordnet ist. wobei S| 
Schwachkontraste angehoben werden, bevor ein Aus- 
gangsbild aus den bei der Zerlegung entstehenden und 
danach modiflzierten Bildern zusammengesetzt wird. 
Die Modifikation der Detailbilder erfolgt dabel sowohl in 
Abhdngigkeit von dem jeweiligen Biteiwert als auch in Jm. 
Abhdngigkeit von dem mittleren Bildwert und/oder der 
Variation der Bildwerte in einem Fenster um das jewei- 
lige Pixel herum. 
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Beschreibung 

[0001] Die Erfindung betrifft ein Verfahren zur Ver- 
arbeitung eines Eingangsbildes mitfblgenden Schritten: 

5 

Zerlegung des Eingangsbildes in eine Folge von 
Detaiibitdern, die jeweils nur einen Teil der Ortsfre- 
quenzen des Eingangsbildes enthalten. 
Modifikation mindestens eines Teils der zu den 
Pixeln wenigstens eines der Oetailbilder gehOren- io 
den Blldwerte, 

Rekonstruktion eines Ausgangsblldes aus den ggf. 
gefilterten Detaiibildern. 

[0002] Das Eingangsbild kann dabei eine Rdntgen- is 
aufnahme sein. doch kfinnen auch andere bei der medi- 
zinischen Diagnostik erzeugte und in elektronischer 
Form vorliegende Bilder verarbeitet werden. wobei 
jedem Pixel (als Pixel wird ein Bildelement bezeichnet; 
ein BiM setzt sich aus einer Matrix von Pixeln zusam- 20 
men) ein die Helllgkeit in diesem Pixel tennzeichnender 
Blldwert zugeordnet 1st. AuBerdem bezleht sIch die 
Erfindung auf eine Anordnung zur Durchfuhrung dieses 
Verfahrens. 

[0003] Ein Verfahren der eingangs genannten Art 25 
ist aus der EP-A 527 525 bekannt. Bei dem bekannten 
Verfahren wird in einer ersten Zerlegungsebene das 
Eingangsbild einer TiefpaOfilterung unterzogen. und 
das entstehende TIefpaBbild wird von dem Eingangs- 
bild subtrahlert, wobel sich ein als erstes Detailbild 30 
bezeichnetes Differenzbild erglbt. In dem nur noch die 
allerhOchsten Ortsfrequenzen des Eirrgangsbikies 
erhalten sind. Wegen der Subtraktlon des Tiefpa3bildes 
von dem Eingangsbild kCnnen sich fOr Pixel in dem 
Detailbild sowohl positive als auch negative Werte erge- 35 
ben. Diese stellen den Kontrast des Bildes dar; doch 
werden sie im folgenden gleichwoN als Bildwerte 
bezeichnet 

[0004] In einer zwelten Zerlegungsebene wind das 
TiefpaBblkj einer erneuten TiefpaOfilterung unterzogen, 40 
so daB ein zwertes TiefpaBbild entsteht. das bezogen 
aufdas Orginalblld weniger hoherfrequente Ortsfre- 
quenzanteile reprdsentiert. als das erste TIefpaBbild. 
Das zweite TiefpaBbild wird vom ersten subtrahiert, 
wobei ein zweites Detailbild entsteht, das QbenA/iegend 45 
Ortsfrequenzen reprdsentiert. die unterhalb des Orts- 
frequenzbandes des ersten DetailbiUes. aber oberhalb 
des Ortsfrequenzberelchs des zwelten TiefpaBblMes 
liegen. 

[0005] In einer drrtten Zerlegungsebene wird auf so 
analoge Weise aus dem zwelten TiefpaBbikd ein drittes 
TiefpaBbild mit noch weniger hOherfrequenten Ortsfre- 
quenzantellen und ein aus der Differenz des zweiten 
und dritten TiefpaBbildes abgeleitetes drittes Detailbild 
erzeugt. Das von diesem dritten Detailbild reprasen- ss 
tierte Ortsfrequenzband liegt Qbenviegend unterhalb 
des Ortsfrequenzberelchs des zweiten Detailbildes. 
aber oberhalb des Ortsfrequenzbereichs des dritten 



TiefpaBbildes. An diese dritte Zerlegungsebene schlie- 
Ben sich noch weitere Zerlegungsebenen an, wobei die 
entstehenden Derailbikler stets niedrlgere Ortsfrequen- 
zen reprasentleren. als das Detailbild der vorhergehen- 
den Zerlegungsebene. 

[0006] Die Detailbilder und das Restbiid werden 
addiert, nachdem mindestens eines der Detailbilder 
modifiziert worden ist. Damit sollen die diagnostisch 
relevanten Bikiinformationen hervorgehoben und die fur 
die Diagnose unwichtigen oder stOrenden BlUlnfbrma- 
tionen unterdrQckt werden. 

[0007] In der EP-A 527 525 sind zwei Altematlven 
for die Modifikation der Detailbilder angegeben: 

a) Bei einer ersten Alternative werden die Biklwerte 
eines Detailbildes Pixel fur Pixel mit einem Faktor 
multipliziert, der von einem Bildwert eines Pixels 
mit entsprechenden Koordinaten in einem TiefpaB- 
bild abhdngt, dessen Ortsfrequenzbereich unter- 
halb des Ortsfrequenzbandes des Detailbildes 
liegt Bel niedrlgen Bildwerten In diesem TiefpaB- 
blki (die bei einer Rdntgenaufnahme den hellen 
BiMpartien entsprechen) Ist der einem Tabellen- 
speicher (look up table) entnommene Faktor kl einer 
als 1 und bei groBen Mittelwerten grdBer als 1. 
Dadurch werden die Kbntraste in den hellen Bildbe- 
reichen. wo sich das Rauschen besonders stOrend 
bemerkbar macht, abgesenkt. wdhrend sie in den 
dunkteren Bildberelchen angehoben werden. Aller- 
dlngs erfolgt die Absenkung bzw. die Anhebung 
unabhdnglg von den BIMwerten Im Detailbild. Das 
Rauschen In diesen BiMpartien wird dadurch also 
in gleicher Weise abgesenkt, wie die in den grOBe- 
ren Bildwerten enthaltene Nutzinfbrmation. 

b) Bei der zweiten Alternative werden die Bildwerte 
wenigstens eines Detailbildes in Abhdngigkeit von 
ihrer GrOBe modifiziert, entsprechend einer nicht 
linearen Kennlinle, wodurch die ganz niedrigen 
Bildwerte (entsprechend dem Rauschen) abge- 
schwacht, mdBIg niedrige Bildwerte angehoben 
und hohere Bildwerte unverdndert bleiben. Aller- 
dings erfolgt diese Modifikation unabhdngig davon. 
ob der Bildwert einem Pixel zugeordnet ist, das in 
den hellen oder in den dunklen Bikfpartien liegt. 
Auch in den dunkleren Bildpartien, in denen sich 
Rauschen in geringerem MaBe auswirkt als in den 
helien Bikipartien. werden also niedrige Bildwerte 
des Detailbildes im Verhaitnis zu hOheren Bildwer- 
ten geschwdcht. Dieses Verfahren hat die Nach- 
telle, daB es zum einen nIcht robust gegenuber 
Rauschpulsen ist, die nicht mehr im Bereich der 
ganz schwachen Kbntraste sondern im Bereich 
mdBig schwacher Kontraste liegen und deshalb 
soger verstdrkt werden und zum anderen Nutzinfor- 
mation im Bereich ganz schwacher Kbntraste in 
jedem Fall geschwdcht wird. 
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[0008] Aufgabe der vorliegenden Erfindung ist es. 
ein Verfahren der eingangs genarmten Art so zu gestal- 
ten. da6 die Bildqualitat noch welter verbessert werden 
kann. Diese Aufgabe wird erfindungsgemaB dadurch 
gelOst, daf3 die Modifikation eines Biidwertes der Muttl- 
plikation dieses Biidwertes mit einem Faktor entspriciit. 
der sowohl von diesem Bildwen selbst als auch von 
dem mittieren Bildwert und/oder der Variation der Bild- 
werte in einem Fenster urn das jeweilige Pixel In dem 
Eingangsbild Oder einem daraus abgeleiteten Bild 
abhangt. 

[0009] Die Erfindung basiert auf der Erkenntnis, 
da3 sich eine optimale Bildqualitdt nur dann ergeben 
kann. wenn die Modifikation der Bildwerte sowohl von 
dem jeweiligen Bildwert selbst als auch von dem mittie- 
ren Bitdwert und/oder der Variation der Bildwerte in 
einem Fenster um das jeweilige Pixel in dem Eingangs- 
bild Oder einem daraus abgeleiteten Bild abhangt 
Dadurch Ist es belspielswelse moglich. da3 in einem 
Detailblld die niedrigen Bildwerte im Verhdltnis zu den 
hflheren Bildwerten dann verstdrkt werden, wenn sie 
einem Pixel zugeordnet sind. das in einer dunklen Bild- 
partie liegt, und nicht verstdrkt werden, wenn das 
betreffende Pixel in einer hellen Bildpartie liegt. 
[0010] Die Erfindung beruht weiterhin auf der 
Erkenntnis, da 3 die Modifikation statt vom mittieren 
BiMwert (was gleichbedeutend ist mit der Dichte des 
Blldes) In einem Fenster um das betreffende Pixel 
herum auch von der Variation der BiUwerte in einem 
solchen Fenster abhdngen kann - jedenfalls dann, wenn 
es darum geht das Rauschen bei der Kbntrastanhebung 
nicht sichtt>ar mit anzuheben. In Bildbereichen, in 
denen die Bildwerte stark variieren. wird das Rauschen 
ndmlich nicht so stark wcUirgenommen, wie in den Bild- 
bereichen. in denen nur eIne schwache Variation der 
Bikiwerte gegeben ist. 

[0011] Eine MOglichkeit, die Erfindung zu imple- 
mentieren bestOnde darin, einen Tabellenspeicher vor- 
zusehen, in dem for jeden Biklwert eines Pixels und fOr 
jede Dichte (mittlerer Bildwert) in einem Fenster um die- 
ses Pixel ein Faktor gespeichert ist. Die Speicherung 
einer solchen zweifachen Abhangigkeit in einem Tabel- 
lenspeicher ware aber aufwendig. und Anderungen die- 
ser zweifachen Abhdngigkeit wdren nur mit 
auBerordentlichem Aufwand zu implementieren. 
[0012] Anspruch 2 beschreibt demgegenuber eine 
wesentlich einlachere MOglichkeit. Dabel werden 
zunachst die niedrigen Bildwerte in dem Detailbild (d.h. 
die schwachen Kbntraste) angehoben. Jedoch wind das 
AusmaB dieser Anhebung in Abhangigkeit von der 
Dichte (und/oder der Variation der Bildwerte) gesteuert. 
Diese Steuerung kann dann so definiert sein, daf3 in 
hellen Bildpartien keine Korttrastanhebung erfolgt, wohl 
aber in dunklen Bildpartien. Dadurch wird das Rau- 
schen zwar nicht verringert; es wird durch die Kbntrast- 
anhebung jedoch auch nicht verstarkt. 
[0013] Eine besonders einfache Art der Steuerung 
ist in Anspruch 3 angegeben. Es erfolgt eine gewichtete 



Summierung des Detailbildes vor und nach der Konver- 
tierung. wobei die Gewichtung von der Dichte (dem 
mittieren Biklwert) und/oder der Variation der Bildwerte 
in dem Fenster um das Pixel im Orglnalbild oder einem 

5 abgeleiteten Bild abhangt. Wenn das Detailbild vor der 
Kbnvertierung mit einem Gewichtungsfaktor 1 und das 
Detailbild nach der Kbnvertierung mit einem Gewich- 
tungsfaktor 0 gewichtet wird, enthait das Summenbild 
keine Kbntrastanhebung. Wird hingegen das Detailbild 

10 vor der Konvertierung mit einem Gewichtungsfaktor 0 
und das Detailbild nach der Kbnvertierung mit einem 
Gewichtungsfaktor 1 beaufschlagt. dann ergibt sich die 
voile Kbntrastanhebung. Aquivalent dazu ist es, (Pixel 
fOr Pixel) die Bildwerte des Detailbildes vor und nach 

75 der Konvertierung voneinander zu subtrahieren und die 
Differenz gewichtet zu dem nicht konvertierten Detail- 
bild zu addieren. 

[0014] Anspruch 4 beschreibt eine fur die Unter- 
drOckung des Rauschens geeignete Weiterbildung, die 

20 davon ai^geht, da3 der wesentllche Anteil des Rau- 
schens sich auf das erste Detailbild mit den hOchsken 
Ortsf requenzen und ggf. noch auf das zweite DetailbiU 
konzentriert. Wenn mit abnehmender Dichte und 
ahnehmender Variation der Bildwerte die Schwachkon- 

25 trastanhebung reduziert wird. wird das Rauschen in den 
dafur empfindlichen Biklbereichen (helle Bildbereiche 
bzw. Bildbereiche mit niedriger Variation der BikJwerte 
bzw. hohe Ortsfrequenzen) nicht angehoben. 
[0015] Hierbei hdngt das AusmaB der Bildverbes- 

30 serung also von der Wahrscheinltehkeit ab, das Bildsi- 
gnal und nicht das Rauschen zu verstarken, bzw. von 
der Wahrscheinlichkelt. daB eventuell mirverstarktes 
Rauschen maskiert wird. Der Vorteil dieser Art der 
rauschrobusten Bildverbesserung ist, daB keine 

35 Abschwachung von Originalbildinhalten auftritt. 

[0016] Anspruch 5 beschreibt eine Vorrichtung zur 
Durchfuhrung des erfindungsgemaBen Verfahrens. Die 
weitere Ausgestaltung nach Anspruch 6 bezleht sich 
auf ein ROntgensystem mit einer Speicheranordnung 

40 (Datenbank) in der for eine Anzahl von Organen je ein 
Satz von Aufnahmeparametern gespeichert ist. ROnt- 
gensysteme mit solchen Datent>anken sind unter dem 
Begriff "APR" (Anatomically Programmed Radiography) 
an sich seit langem bekannt. Dabei wird von der 

45 Erkenntnis ausgegangen, daB es keine Einstellung der 
Aufnahmeparameter gibt. die fur die Aufnahmen sdmtii- 
cher Organe optimal ist. Da bei einer Rdntgenaufnahme 
dieser Art die Aufnahmeparameter ohnehin organab- 
h&igig bzw. In Abhdngigkeit von der darzustellenden 

so anatomischen Region des Patienten aufgenommen 
werden, sieht diese Weiterbildung vor, daB auBer den 
erforderlichen Aufnahmeparametern auch die fOr die 
ROntgenaufnahme der betreffenden anatomischen 
Region optimalen Bildverarbeitungsparameter gespei- 

55 chert werden. 

[001 7] Die Erfindung wird nachstehend anhand der 
Zeichnungen naher eriautert. Es zeigen: 
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Fig. 1 ein BOntgenaufnahmegerdt mit einer 

erfindungsgemdBen Bildverarbeitungs- 
vorrichtung. 

Fig. 2 ein Blockdiagramm zur Erlduterung der 

Bildverarbeitung, 
Fig. 3 ein Blockdiagramm der darin enthalte- 

nen Modrfikationsstufe, 
Fig. 4 eine Kennlinle zur Anhebung der schwa- 

Chen Kontraste. 
Fig. 5 a 5b verschiedene Kennlinien zur Steuerung 

der Schwachkontrastanhebung. 
Fig. 6 ein Dichte/Ortsfrequenzdiagramm und 

Fig. 7 eine weitere Ausgestaltung der Modifi- 

kationsstufe. 

[0018] In Fig. 1 ist mit 1 ein ROntgenstrahler 
bezeichnet, der ein Strahlenbundel 2 emittiert. das 
einen nicht ndher dargesteilten, auf einem durch eine 
Tischplatte 3 angedeuteten Parientenlagerungstisch 
gelagerten Patienten durchstraNt und auf einen Bildauf- 
nehmer 4 trifft, der die auftreffende Rdntgenstrahlungs- 
intensitdt onsabhdngig in elektrische Signale umsetzt. 
[0019] Der BiMaufnehmer 4 kann beispielsweise 
aus einer Vielzahl von matrixfOrmig angeordneten. fOr 
ROntgenstrahlung empfindlichen Detekrorelementen 
bestehen. die elektrische Signale erzeugen. deren 
Amplitude der Inlensitat der auf sie auftreffenden ROnt- 
genstrahlung entsprechen. Der Bildaufnehmer kann 
aber auch einen Photoleiter enthalten, auf dem ein von 
der ROntgenstrahlung abhangiges Ladungsmuster 
erzeugt wird, das elektrosratisch abgetastet wird. 
SchlieGlich kann als Bildaufnehmer auch ein Speicher- 
phosphor vorgesehen sein, in dem durch die ROntgen- 
strahlung ein latentes Bitd erzeugt wird, das mittels 
eines Lasers ausgelesen und digiralisiert werden kann. 
Die konkrete physikalische Ausgestaltung des Bildauf- 
nehmers ist fur die Erfindung nicht wesentiich; wichtig 
ist nur, da6 das ROntgenstrahlenbild in eine Folge von 
digitalen Datenworten umgesetzt werden kann. die die 
Bildwerte an den einzelnen Pixein eines Eingangsbildes 
reprdsentieren. 

[0020] Der Bildaufnehmer 4 wird von einem daran 
angepaf^en Lesegerm 5 ausgelesen, das eine Folge 
von Digitalwerten erzeugt, die jeweils einem Bildwert 
des Eingangsbildes entsprechen. Das Lesegerm 5 
kann dabei bereits so ausgestaltet sein, da6 es die 
durch den Bildaufnehmer bedingten Artefakte beseitigt 
und beispielsweise schon eine Korttrast-Harmonisie- 
rung vornimmt, durch die der Dynamikumfang der gro- 
ben Strukhjren in diesem Bild reduziert wird. 
[0021] Die von der Leseeinhert 5 erzeugten Digital- 
werte werden einer Workstation 6 zugefOhrt. in der 
einerseits eine Bildverarbeitung erfdgt andererseits 
aber auch eine Steuerung eines ROntgengenerator 7, 
an den der ROntgenstrahler 1 angeschlossen ist. Die 
Workstation 6 arbeitet mit einem Monitor 8 zusammen, 
auf dem ein ROntgenbiM wiedergegeben werden kann, 
das aber auch den Systemzustand wiedergeben kann. 



z.B. die im Rdntgenstrahler wirksame ROhrenspan- 
nung, die Aufnahmedauer oder die KOrpen-egion, die 
bei der ndchsten ROntgenaufnahme abgebildet werden 
soil. Au8eidem ist eine Eingabeeinheit 9 vorgesehen 

5 mit der der Benutzer Steuerbefehle eingeben kann. 
SchlieBlich wirkt die Wbrkstation 6 mit einer in einem 
Speicher 10 enthattenen APR-Datenbank zusammen. 
die einerseits die Aufnahmeparameter und andererseits 
Steuerparameter for eine nachfolgende Bildverarbei- 

10 tung jeweils fur verschiedene Organe bzw. verschie- 
dene anatomische Regionen enthalten kann. 
[0022] Nachfblgend wird die in der Workstation 6 
durchgefOhrte Verarbeitung eines Eingangsbildes Iq 
anhand des Blockdiagramms gemdB Fig. 2 ndher erldu- 

15 ten. Das Eingangsbild Iq wird zundchst im Block 61 1 
einer TiefpaBfilterung unterzogen, wobei die oberste 
Ortsfrequenz die nach der TiefpaBfilterung ubrig bleibt 
halb so gro3 ist wie die maximale Ortsfrequenz, die in 
dem Eingangsbild Iq enthalten sein kann. Gleichzeitig 

20 erfblgt, wie durch den abwdrtsweisenden Pfeil ange- 
deutet, eine Unterabrastung. indem das TiefpalSbild in 
Zeilen- und in Spaltenrichtung nur for jeden zweiten 
Bildpunkt berechnet wird. Das auf diese Weise 
erzeugte TiefpaBbild ^ wind u.a. einer Stufe 621 zuge- 

25 fuhrt. in der einerseits die Zahl der Bildpunkte vergrO- 
Bert wird, in dem beispielsweise zwischen je zwei 
Zeilen bzw. Spalten des TiefpaBbiides 1^ eine Zeile bzw. 
eine Spalte von Nullen eingefQgt wird und das entste- 
hende Bild mit einem TiefpaBkernet gefaltet wird. Das 

30 Bild am Ausgang der Stufe 621 entspricht daher im 
wesentlichen dem Bild. das sich ergeben wQrde, wenn 
in der Stufe 61 1 nach der TiefpaBfilterung nicht jede 
zweite Zeile und Spalte eliminiert worden ware. Jedoch 
hat dieses TiefjpaBbild genauso viele Pixel wie das Ein- 

35 gangsbild Iq, aus dem es abgeleitet ist. Im Block 631 
wird dieses TiefpaBbild - Pixel fur Pixel - von dem Ein- 
gangsbild Iq subtrahiert. Das sich dabei ergebende 
Subtraktionsbild stellt das erste Detailbitd S^ dar, das 
lediglich die hoheren Ortsfrequenzen des Bildes Iq ent- 

40 halt, aus dem es aksgeleitet wurde. 

[0023] In einer zweiten Zerlegungsebene wird das 
TiefpaBbild 1^ , welches die urrtere Haifte des Frequenz- 
spektrums des Orginalbiides reprasentiert. in der glei- 
chen Weise verarbeitet, wie in der vorangehenden 

45 Zerlegungsebene das Eingangsbild Iq. Das Bild wird 
also in einer Stufe 612 einer TiefpaBfilterung unterzo- 
gen, und gleichzeitig wird die Zahl der Pixel in Zeilen- 
und in Spaltenrichtung wiederum halbiert. so daB sich 
ein TiefpaBbild I2 mit erneut um einen Faktor 2 reduzier- 

so ter Zeilen- und Spaltenzahl erglbt. welches bezogen auf 
das Orginalbild nur noch das untere Viertel des Fre- 
quenzspektrums enthait. Das TiefpaBbild I2 wird in 
einer Stufe 622 wiederum so interpoliert, daB sich ein 
TiefpaBbild mit der gleichen Zahl von Pixein wie das 

55 TiefpaBbild l^ ergibt. Dieses hinter dem Block 622 
erzeugte TiefpaBbild wifd in dem Block 632 von dem 
Bild I1 subtrahiert. aus dem es abgeleitet wurde, so daB 
sich ein zweites Detailbild S2 ergibt. Dieses Detailbitd 
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reprasentiert bezogen auf das Orgnialbild die obere 
Hdlfte des Ortsfrequenzbandes, das sich an das im 
Detailbild enthattene Ortsfrequenzband anschliefSt. 
umfaet aber nicht den im TjelpaBbild I2 enthaltenen 
Ortsfrequenzbereich. 

[0024] In einer dritten Zerlegungsstufe wird das 
Tiefpa3bitd I2 in derselben Weise verarbeitet, wie in der 
zwetten Zeriegungebene das Tiefra3bild 1^ Es resultie- 
ren daraus ein Tiefpa(3bild I3 (mit der haiben Zeilen- und 
Spattenzaht wie das TiefjpaBbild I2}. welches das untere 
Achtel des Orginalspektrums reprdsentiert sowie ein 
Detailbild S3, dessen Ortfrequenzband sich an das 
Ortsfrequenzband des Detallbildes S2 anschlieBt. Der 
ZerlegungsprozeB kann n solcher identisch aufgebau- 
ter Zerlegungsebenen enthalten. Typische Werte fOr n 
sind 7 Oder 8. Das bedeutet. daB das TiefpaBbild Ip (das 
auch ats Restbild bezeichnet werden kann) 128 oder 
256 Zeilen bzw. Spatten im Eingangsbild durch eine ein- 
zige Zeiie bzw. Spatte ersetzt. 

[0025] BeKanntlich kOnnte man das Eingangsbild Iq 
vollstandig aus dem Restbild und den Detailbildern 
Si...Sn rekonstruieren, wenn man die Bilder zueinander 
addiert, nachdem man durch Interpolation dafOr 
gesorgt hat. daB die zu addierenden Bilder gleiche Zei- 
len- und Spaltenzaht aufweisen. Diesem Zweck dienen 
die in den Zerlegungsebenen enthaltenen Interpotati- 

onsstufen 641 64n und die ihnen in den einzetnen 

Ebenen nachgeschalteten Additionsstufen 651, 652. 
653 und 65n. Dieser RekonstruktionsprozeB ist in dem 
eingangs enwdhnten Dokument im einzelnen beschrie- 
ben. worauf zur Vermeidung von Wiederholungen 
Bezug genommen wird. 

[0026] Da aber mindestens eines der Detaiibilder 
Si....Sn durch eine Modifikationsstufe M-i. Mg, M3....Mn 
modifiziert wird, so daB modifizierte Detaiibilder Smi. 
Sm2.Sm3 -Sn entstehen, weicht das durch die Rekon- 
struktion entstehende Ausgangsbild Imvon dem Ein- 
gangsbild Iq ab. In Fig. 3 ist eine Modifikationsstufe Mj 
(i= 1,2.3 n) dargestellt. 

[0027] Der Biklwert eines jeden Pixels im Detailbild 
S; wird in einer Stufe 21 mit nichtlinearer Gbertragungs- 
kennlinie transformiert. Dies erfolgt am einfachsten mit- 
tels eines label lenspeichers (look up table) . der jedem 
Bildwert Sj an seinem Eingang einen ggf. gednderten 
Bildwert Sgj an seinem Ausgang zuordnet. Fig. 4 stellt 
den typischen Verlauf der nichtlinearen Kennlinie (mit 
dem Bildwert Sj als Abszisse und dem Bildwert Sei als 
Ordinate) gemflB Block 21 dar, doch k6nnen die Breite 
des angehobenen Bereiches und das AusmaB der 
Anhebung bei den verschiedenen Modifikationsstufen 
unterschiedlich groB sein. Man erkennt, daB groBe Bild- 
werte unverdndert bleiben. wdhrend Meinere Biklwerte 
(die im ersten Detailbild dem Rauschen entsprechen 
kOnnen) relativ zu den gr6Beren Bildwerten angehoben 
werden. Da die Bitdwerte der Detaiibilder den Kbntrast 
des Eingangsbildes in dem jeweiligen Ortsfrequenzbe- 
reich darstellen. bewirktdie angegebene Ubertragungs- 
kennlinie eine Anhebung der schwachen Kontraste In 



dem jeweiligen Bild. 

[0028] AnschlleBend wird der Bildwert eines jeden 
Pixels Im Detailbild Sj in einer MuKiplikationsstufe 20 mit 
einem Faktor multiplizlert der grOBer ist als 1. Das 

5 bedeutet. daB dieses Detailbild Sj im Verhditnis zu dem 

Restbild In (wenn die Modifikationsstufen Mj^i Mn 

diese Multiplikation nicht enthalten aber auch im Ver- 
hdJtnis zu den Detailbildern der nachfolgenden Zerle- 
gungsstufen Sj^.i...Sn) angeholjen werden. Diese 

10 Multiplikation und die damit einhergehende Betonung 
der hOheren Ortsfrequenzen kann aber ggf. auch entfeil- 
len Oder bereits in der Leseeinheit 5 vorgenommen wer- 
den. 

[0029] Danach wird in der Stufe 22 die Differenz A 

15 zwischen dem Bildwert am Ausgang der Stufe 21 und 
dem im gleichen Pixel in dem Detailbild Sj zugeordne- 
ten Bildwert gebildet. Diese Differenz A wird in einer 
Multiplizierstufe 23 mit einem Faktor b multipliziert, der 
zwischen 0 und 1 liegt. Der multipliziert e Differenzwert 

20 bA wird zu dem Bildwert Sj addiert (Block 24) und bildet 
damit einen Bildwert des modif izierten Detallbildes S^j. 
Wenn der Faktor b^o ist. entspricht der BiMwert nach 
der Addition im Block 24 genau dem Wert im Detailbild 
Sj. Die Kontrastanhebung im Block 21 kommt dabei 

25 nicht zur Wirkung. Wenn b hingegen den Wert 1 hat. ist 
die Kontrastanhebung voll wirksam. 
[0030] Der Faktor b wird aus einem Bild abgeieitet, 
das zu jedem Pixel im Detailbiki Sj einen der mittleren 
Helligkeit an dieser Stelle entsprechenden Wert enthait. 

30 Zum Beispiel kann hierfOr das Eingangsbild Ij^i venA^en- 
det weiden, aus dem das Detailbik:! Sj abgeieitet ist, weil 
dieses BiM genau soviete Pixel enthait wie das Detail- 
bild. Dieses Bild wird einer TiefpaBfilterung unterzogen. 
wobei der mittlere Bildwert - d.h. die Dichte D - in einem 

35 Fenster von z.B. 5 x 5 Pixein herangezogen wird, in des- 
sen Mitte sich das gleiche Pixel befindet, dessen Bild- 
wert jeweils der Schwachkontrastanhebung im Block 21 
unterzogen wird. Alternativ kann aber auch - wie gestri- 
chelt angedeutet - Bild Ij' aus dem Rekonstruktionspfad 

40 (d.h. vom Ausgang der Interpolationsstufe 64j ..i^i .2,3) 
zur Steuerung venwendet werden, wobei die TiefpaBfil- 
terung in Block 25 entfollen kann. Dies ist insbesondere 
dann zweckmdBig, werm das Restbild In durch eine 
zusdtzliche Dichtekorrektur-Stufe 660 einer Dichtekor- 

45 rektur (z.B. entsprechend einer Filmkurve) unterzogen 
wird. Der der Dichte D ensprecherxje gefilterte Bildwert 
wird einem Tabellenspeicher LUTj zugefQhrt der jedem 
Wert D einen Wert b zuordnet. 
[0031] Wie stark die Kontrastanhebung wirksam 

so wird, hangt somit von der Dichte in einem Fenster urn 
das jeweilige Pixel ab sowie von dem Verlauf der in dem 
Tabellenspeicher LUTj gespeicherten Kennlinie. Im fol- 
genden wird beschrieben, wie man die Rauschverstdr- 
kung beschrdnken, aber auch daruberhinaus die 

55 Kontrastanhebung organspezifisch steuern kann. Letz- 
teres funktioniert unter der Voraussetzung. daB eine 
Unterscheidung zwischen Organen mit Hilfe der Krite- 
rien Dichte und DetailgrOBe mOglich ist. Eine solche 
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organspezifische Steuerung der Kontrastanhebung ist 
insbesondere bei Lungenaufnahmen mOglich und wun- 
schenswert. 

[0032] Zur nSheren Erlduterung wird zundchst auf 
Fig. 6 verwiesen. die ein Diagramm darstellt mit der 
Dichte D als Ordinate und den Ortsfrequenzbdndern 
der einzelnen Detaiibilder als Abszissen. In diesem Dia- 
gramm ist eine Region dargestellt. die sich in den 
Detailbildern , S2 mit den hOchsten Ortsfrequenzen 
bei der niedrigsten Dichte befindet. Der zugehdrige 
Dichtebereich entspricht in einem ROntgenbild einer 
niedrigen Strahlendosis. Wegen des nichtlinearen, im 
aligemeinen logarithmischen Zusammenhanges zwi- 
schen der Strahlendosis und der Dichte. wird bekannt- 
lich in diesem Bereich das Rauschen am stdrksten 
sichtbar. 

[0033] AuBerdem ist in Fig. 6 ein zweiter Bereich R2 
eingezeichnet. der einer mittleren Dichte bzw. einer 
mittleren Helligkeit entspricht und der in einem niedrige- 
ren Ortsfrequenzbereich liegt, z.B. in den Detailbildern 
S3 und S 4. In dieser Region werden bei einer Lungen- 
aufnahme mit heute ubiichen Aufl6sungen von 0,2mm/ 
ein GroBteil der Lungengefd3e abgebildet. Werden in 
diesem Bereich die Schwachkontraste angehoben, 
dann kann sich auch bei gesunden Lungen ein patholo- 
gischer BildeindrucK ergeben. Hingegen sollen im glei- 
chen Ortsfrequenzbereich die Strukturen im unteren 
Dichtebereich in ihrem Kontrast angehoben werden. In 
diesem Bereich liegen z.B. Strukturen. die von Media- 
stinum und Herzschatten Qberlagert werden und daher 
schlecht sichtbar sind. Gleichzeitig spielt in diesem 
Ortsfrequenzbereich das Rauschen keine so gro3e 
Rolle. 

[0034] In Fig. 5a und 5b sind die dafOr erforderli- 
Chen Kennlinien dargestellt, die in den Tabellenspei- 
chern LUTj der Modtfikationsstufen gespeichert sind. 
[0035] Rg. 5a stellt die Steuerkennlinie fur das 
hOchste Ortsfrequenzband S^(und ggf. 82) dar. Bei 
niedriger Dichte ergibt sich dabei ein niedriger Wert des 
Faktors b, d.h. die durch die Kennlinie gemSB Fig. 4 
bewirkte Kontrastanhebung wird am Ausgang der Stufe 
24 kaum wirksam. Die in dem Bereich von R^ liegenden 
Schwachkontraste - insbesondere das Rauschen -wer- 
den also nicht angehoben werden. Das Rauschen wird 
nicht reduziert, praktisch jedoch auch nicht verstdrkt. 
wie es der Fall wdre. wenn alle Schwachkontraste ange- 
hoben wQrden. Zu grOBeren Dichtewerten hin kann der 
Faktor b zunehmen und damit das AusmaB in dem die 
Kontrastanhebung wirksam wird. In diesem Dk:htebe- 
reich ist das Rauschen im aligemeinen nicht so stOrend 
wie im unteren Dichtebereich. 
[0036] Fig. 5 b zeigt demgegenQber eine Kennlinie. 
wie sie in dem Tabellenspeicher LUT3 fur die Modifika- 
tion im Block M3 enthalten sein kann. Bei niedrigen 
Werten der Dichte ist der Faktor b groB. Dadurch wird 
die Schwachkontrasranhebung bei niedriger Dichte im 
ROntgenbild (also im Bereich des Mediastinums, des 
Herzschattens und des Abdomens) voll wirksam. Zu 



h6heren Werten der Dichte hin nimmt der Faktor b 
gemdB Fig. 5b kontinuierlich ab, so daB in dem Dichre- 
bereich, in dem die Lungengef^Be uberwiegend abge- 
bildet werden, die Schwachkontrastanhebung in weit 

5 geringerem MaBe wirksam wird. 

[0037] Die Kennlinie gemdB Fig. 5b ist somit fQr 
Lxingenaufnahmen gut geetgnet. FQr andere Aufnah- 
men hingegen kann ein vOllig anderer Kennlinienverlauf 
erforderlich sein. Die fur verschiedene Organe optima- 

10 len Kurven bzw. Bildverarbeitungsparameter sollten 
daher - zumindest fQr einige ihrer StOtzpunkte - in der 
APR Datenbank 10 gespeichert sein und bei einer Vor- 
gabe des darzustellenden Organs durch den Untersu- 
Cher aufgerufen bei der nachfblgenden Bildverabeltung 

IS vorgegeben werden. 

[0038] Fig. 7 zeigt eine andere Ausfuhrungsform 
der Modifikationsstufe Mj. Das Detailbild Si wird dabei 
wiederum im Block 21 einer Schwachkontrastanhebung 
unterzogen. Das resultierende Detailbild mrt verstdrk- 

20 tem Kontrast Sej wird in einer Kbmbinationsstufe 27 mit 
dem Detailbild Sj BildpunktfQr Bildpunkt kombiniert. FQr 
das modlf Izierte Ausgangsbild gilt dann die Gleichung 

Smi=b-Se,+(1-b)-S, (1) 

25 

[0039] Diese gewichte Summierung von S^ und Sqi 
ist Equivalent zu der aus Fig. 3 ersichtlichen Kombina- 
tion, wobei zundchst die Differenz A zwischen dem im 
Kontrast verstdrkten Detailbild Sei und dem Detailbild Sj 

30 selbst geblMet wird und die Differenz A dann mit dem 
Qewicht b zu dem Detailbild Sj addiert wird. 
[0040] Der Wert b ergibt sich dabei als Produkt der 
Faktoren b^ und b^. wobei b^ auf die gleiche Weise 
ermittelt wird wie der Faktor b gemdB Fig. 3. Der Faktor 

35 bs hangt von der Variation der Bildwerte in einem Fen- 
ster um das Pixel ab, das gerade gemdB Gleichung 1 
bearbeitet wird. Bei einer schwachen Variation der Bild- 
werte ist der Faktor b^ klein, und bei einer starken Varia- 
tion ist der Faktor bg groB. Damit wird der Tatsache 

40 Rechnung getragen, daB das Rauschen bei einem 
Pixel, in dessen Umgebung sich die BiMwerte stark 
dndern, weniger stOrend ist als bei einer Pixel, in des- 
sen Umgebung die Bildwerte nur geringfugig schwan- 
ken. Die Variation der Bildwerte kann dabei auf 

45 unterschiedlichste Weise bestimnrt werden. Beispiel- 
haft seien hier drei verschiedene VariationsmaBe ange- 
geben. 

[0041] Als erstes Beispiel sei hier die Standardab- 
weichung 

so 



T= jj2:(<BW(Sj))-BW(Si))' (2) 

55 

angegeben. Hierbei bezeichne BW(Si) einen Bildwert 
von Sj und <BW(S|)> den Mittelwert der Bildwerte in 
einem lokalen Fenster um einen betrachteten Pixel. C 
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steht fur aile Pixel, die sich in einem tokalen Fenster der 
Gr6Be F um das jeweilige Pixel befinden. Neben der 
Standardabweichung sind natQrIich auch alie anderen 
Indikatoren der statistischen Schwankungsbreite, z.B. 
die hOheren Momenta, geeignete Ma6e f Or die Variation 5 
der Bildwerte. QleichermaBen geeignet sind auch histo- 
gram-basierte MaBe die nach folgender Beziehung 
bestimmt werden 

T=-2:*(Pg) (3) 

9u 



[0042] Dabei ist Pg der Bruchteil der Pixel mit dem 
Bildwert g, und gu und go stellen den minimalen bzw. 
maximalen Grauwert dar. f(pg) ist eine strikt konvexe 
Funktlon. fur die z.B. die Beziehung gelten kann 

Oder die Beziehung 

f(Pg) = Pg'' (5) 

[0043] Bei den beispielhaft genannten Variations- 
maBen T (Gleichungen 2 und 3) fuhrt eine geringe 
Variation der Bildwerte zu einem kleinem Wert fOr T und 
eine groBen Variation zu einem groBen Wert. 
[0044] Nach der Ermittlung des Wertes T fur das 
jeweilige Pixel gemdB den Gleichungen 2 bzw. 3 und 4 
Oder 5 im Block 28 wird dem jeweils ermittelten Wert T 
mit einem Tabellenspeicher 29 ein Wen bg zugeordnet 
der groB ist fOr groBe T und klein fOr kleine Werte von T 
Die Werte fur bj und bs werden dann im Block 30 mit- 
einander multipliziert, woraus sich der Faktor b fur das 
jeweilige Pixel ergibt. Der wesentliche Urrterschied 
gegenQber Fig. 3 besteht darin, daB auch der Wert bs in 
den Faktor b eingeht. 

[0045] Es ist aber auch mOglich, den Faktor b aus- 
schlieBlich von der Variation der Bildwerte in einem 
Fenster um das jeweilige Pixel herum abhdngig zu 
machen. Die BlOcke 25 und 26 kOnnten dann entlallen. 
Die Unempfindlichkeit gegenQber Rauschen wdre dann 
auch bei anderen medizinischen Bildern (z.B. Magnet- 
resonanz-Bildern) gegeben. unabhangig davon. ob 
dabei das Rauschen in hellen Bildpartien starker sicht- 
bar Ist als In dunklen. 

PatentansprOche 

1. Verlahren zur Verarbeitung eines Eingangsbikles 
mitfolgenden Schritten: 

Zerlegung des Eingangsbitdes in eine Folge 
von Detailbildern, die jeweils nur einen Teil der 
Ortsfrequenzen des Eingangsbildes enthalten. 
Modif ikation mindestens eines Teils der zu den 



Pixein wenigstens eines der DetailbikJer gehO- 
renden Bildwerte. 

Rekonstruktion eines Ausgangsbildes aus den 
ggf. gefilterten Derailbildern, 
dadurch oekennzeichnet. daB die Modtfikation 
eines Bildwertes der Multiplikation dieses Bild- 
wertes mit einem Faktor entspricht, der sowohl 
von diesem Bildwert seit>st als auch von dem 
mittleren Bildwert und/oder der Variation der 
Bildwerte in einem Fenster um das jeweilige 
Pixel in dem Eingangsbild Oder einem daraus 
abgeleiteten Bild abhftngt. 

2. Verlahren nach Anspruch 1 

dadurch oekennzeichnet. daB die Modifikation fbl- 
gende Schritte umfaBt: 

Konvertierung der Bildwerte in wenigstens 
einem der Detailbilder mit einer nichtlinearen 
Kennlinie derart. daB die niedrigen Bildwerte in 
dem Detailblld im Vergleich zu den grOBeren 
angehoben werden, 

Steuerung des AusmaBes der Anhebung der 
niedrigen Bildwerte in den PIxeIn des Detailbil- 
des in Abhangigkeit von dem mittleren Bildwert 
und/oder der Variation der Bildwerte in einem 
Fenster um die jeweiligen Pixel in dem Ein- 
gangsbite:! Oder einem daraus abgeleiteten Bild. 

3. Verlahren nach Anspruch 2, 

dadurch oekennzeichnet. daB die Steuerung des 
AusmaBes der Anhebung des Bildwertes eines 
jeden Pixels die gewlchtete Summierung der Bild- 
wertes in dem Detailbild und in dem konvertierten 
Detailbild mit einem Anteilsfaktor umfaBt, der von 
dem mittleren Bildwert und/oder der Variation der 
Bildwerte in einem Fenster um das jeweilige Pixel in 
dem Eingangsbild Oder einem daraus abgeleiteten 
Bild abhangt. 

4. Verfahren nach Anspruch 3, Insbesondere fOr eine 
ROntgenaufnahme als Eingangsbild, dadurch 
gekennzeichnet. daB zumindest fOr die Detailbilder 
mit den hOchsten Ortsfrequenzen die Antellsfakto- 
ren solche Werte haben. daB das AusmaS der 
Anhebung Anhebungen mit at>nehmender Dichte 
und/oder abnehmender Variation reduziert wird. 

5. Vorrlchtung zur Durchfuhrung des Verlahrens nach 
Anspruch 1 mit einem Bildgebungssystem (1, 4. 5) 
und einer Biklverabeitungseinheit (10), die verse- 
hen Ist mit 

- Mitteln (611. 621. 631 61n, 62n. 63n) zur 

Zerlegung des Eingangsbildes (Iq) in eine 

Folge von Detailbildern (Si. S2 SJ, die 

jeweils nur einen Teil der Ortsfrequenzen des 
Eingangsbildes (Iq) enthalten. 
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Mittein zur Modtfikation (M^, M2 Mn) minde- 

stens eines Teils der zu den Pixein wenigstens 
eines der Detailbilder gehOrenden Bildwerte, 

- MHteln (641 . 651 64n, 65n) zur Rekonstruk- 

ton eines Ausgangsbildes (Im) aus den ggf. s 
modrfiziertenen Detailbildern (Smi» Sma 

dadurch qekennzeichnet. da3 die Mittein zur 

Modifikation (Mi, M2 M^) so ausgebildet 

sind, da8 die Modifikation eines Bildwertes der 10 
Muitiplikation dieses Bildwertes mit einem Fak- 
tor entspricht. der sowohl von diesem Bildwert 
selbst als auch von dem mittleren Bildwert 
und/oder der Variation der BiMwerte in einem 
Fenster um das jeweilige Pixel in dem Ein- is 
gangsbild Oder einem daraus abgeleiteten Bild 
abhdngt. 

Vorrichtung nach Anspruch 5. 

dadurch oeke nnzeichnet. daB die Vorrichtung als 20 
Bildgebungssystem einen ROntgenstrahler (1) und 
einen ROntgen-Bildaufnehmer (4) zwecks Darstel- 
lung unterschiedlicher Organe durch ROntgenauf- 
nahmen enthdit. da 3 der ROntgenstrahler von 
einem ROntgengeneraror (7) gesteuert wind und 25 
da6 eine Speicheranordnung (10) vorgesehen ist, 
in der fur eine Anzahl von Organen je ein Satz von 
Aufnahmeparametern gespeichert ist, daB jeder 
Satz neben Aufnahmeparametern for den ROntgen- 
generator Parameter fur die Modifikation der Detail- 30 
bilder in Abhdngikeit von dem mittleren Bildwert 
und/oder der Variation der Bildwerte in einem Fen- 
ster um das jeweilige Pixel in der Eingangsbild oder 
einem daraus abgeleiteten Bild enthdit, und daB bei 
einer Vorgabe eines Organs die fOr dieses Organ as 
gespeicherte die Modifikation der Detailbilder auf- 
gerufen und der Bildverabeitungseinheit vorgege- 
ben wird. 
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Noise reduction in an image 



PCT/IB9S/00850 



The invention relates to a method of processing an image. The invention 
also relates to an image processor and to an x-ray examination apparatus. 

The article Image filtering using multiresolution representation in the 
5 IEEE Transactions of Pattern Analysis and Machine Intelligence, 13(1991)426-440, by S. 
Ranganath, discloses a method for processing an image, wherein the image is decomposed 
into detail images at successive resolution levels, and the detail images are filtered so as to 
generate filtered detail images. The processed image is reconstructed from the filtered detail 
images. 

10 The cited reference discloses image filtering using multiresolution 

representation. Image filtering is implemented as the convolution of the image and a two- 
dimensional mask. Preferably, a pyramid strucmre is applied for the multiresolution 
representation. In the cited reference the detail images of the multiresolution representation 
implemented by the pyramid structure are named 'coded images'. The cited reference 

15 mentions that the *coded unages' are filtered. The known method aims at achieving a good 
compromise between noise removal and loss of resolution. It has been found that the known 
method is not satisfactory with respect to reducing noise while retaining edges in the image. 
In particular, the known method is not suitable for producing a processed medical image 
having a high diagnostic quality, i.e. an image in which small details having low contrast are 

20 displayed so as to properly visible. 

It is an object of the invention to provide a method of processing an 
image, which reduces noise and retains relevant image information more effectively than the 
2S known method. 

This object is achieved by the method of processing an image according to 
the invention, which method comprises the steps of 

decomposing the image into detail images at successive resolution levels, 
filtering detail images to generate filtered detail images. 
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filtering a current detail image in dependence on at least one subsequent detail 
image, said subsequent detail image having a resolution level which is lower than the 
resolution level of the current detail image, and 

reconstructing a processed image from the filtered detail images. 
5 The detail images at successive resolution levels contain image 

information at respective spatial scales. The detail images are ordered according to 
decreasing resolution levels and increasing coarseness, that is the subsequent detail image 
contains coarser image features than the current detail iniage. Fine details are prominent in 
the detail image having the finest (highest) resolution level; progressively, coarser features 
10 become more prominent at the coarser (lower) resolution levels. Preferably, the 
decomposition of the image into detail images is carried out by means of a pyramid 
strucmre. In such a pyramid structure, the finest resolution level is at the base of the 
pyramid, the coarsest resolution level is at the top of the pyramid. In the pyramid strucmre, 
the detail images are derived as frequency-band images, i.e. spatial- frequency band-pass 
15 filtered versions of the image. The filtering involves a local averaging of the detail images so 
as to reduce differences between pixel-values that are due to noise. As an alternative, the 
filtering may include a non-linear transformation of pixel- values. Pixel-values represent 
brightness values such as grey-values or colour- values of the image. Because the filtering of 
the current detail image depends on one or several detail images at coarser resolution levels, 
20 the method according to the invention distinguishes between image structures that appear at 
several resolution levels and image strucmres at only one or a few resolution levels. For 
example, pixel-values of the current detail image are retained when there is a significant 
corresponding structure at coarser resolution levels, but pixel-values of the current detail 
image are locally smoothed when there is hardly any significant strucmre at coarser 
25 resolution levels. In particular, a steep edge in the image appears at almost each of the 

resolution levels, but an extended slope in the pixel- values of the image appears almost only 
at one or a few subsequent resolution levels. It has been found that filtering in dependence on 
several lower resolution levels produces good results when said filtering is controlled by a 
sum of gradients of pixel-values of detail images having several resolution levels. The 
30 method of the invention retains relevant image information, notably directional structure-like 
edges in the image, and reduces noise even in portions of the image containing variations of 
pixel-values on a coarse scale. Because the filtering is adjusted on the basis of image 
information with spatial frequencies lower than the spatial frequencies of the current detail 
image, the filtering is notably insensitive to noise. The method of processing an image is 
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particularly suitable for processing a medical diagnostic image, such as an x-ray image, a 
magnetic resonance image or an ultrasound image. Such a processed medical diagnostic 
image has a very small noise component and renders small details having little contrast 
properly visible. 

5 

Advantageous implementations of a method according to the invention are 
defined in the dependent Claims. 

A preferred embodiment of a method of processing an image according to 
the invention controls the filtering on the basis of the next detail image. Image structures, 
10 such as edges, relating to relevant image information appear most strongly at the next coarser 
resolution level. Since filtering of the current detail image is controlled on the basis of the 
correlation between image information in the current and the next detail images, it is 
effectively achieved to reduce noise while retaining relevant image information in the current 
detail image. 

15 A preferred embodiment of a method of processing an image according to 

the invention controls filtering on the basis of differences between pixel-values in the 
subsequent detail image, image strucmres are predominantly represented by differences 
between pixel- values. Notably image structures in the form of edges gives rise to differences 
between pixel-values at different positions in the image. The current detail image is strongly 

20 filtered, i.e. smoothed, at locations corresponding to locations in a detail image of the 

subsequent resolution level where there is a relatively small local difference between pixel- 
values, especially a small gradient of pixel-values. The pixel-values of the current detail 
image are retained at locations corresponding to locations in a detail image of the subsequent 
resolution level where there is a relatively large local difference between pixel-values. 

25 especially a large gradient of pixel-values. Image structures are predominantly represented by 
oriented structures like edges or lines. These structures have a strong variation in intensity in 
one direction in the image whereas the intensity variation in the direction at right angles 
thereto is small. In the vicinity of such an image structure, the gradient of pixel-values has a 
high magnimde and a direction perpendicular to the orientation of the image structure. Local 

30 averaging perpendicularly to the gradient, i.e. along the elongate image structure hardly 
affects the strength of visually important intensity transitions, notably variations along the 
image structures are mainly due to noise. 

Edges or other steep transitions result in strong local gradients in detail 
images at several resolution levels. Moreover, the direction of maximum local variation and 
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therefore also the direction of the gradient will be substantially equal in detail images at 
several resolution levels. The amplitude of the gradient may vary from level to level 
depending on the steepness of the transition and the size of the structure. Therefore the 
information contained in the gradients for different resolution levels is similar for adjacent 

5 resolution levels. The information in the coarser resolution levels is less affected by noise 
and is therefore more reliable. A strong image structure in the subswquent, preferably the 
next, detail image controls the filter at the current detail image to preserve the pixel- value so 
as to retain relevant image information in the current detail image. 

In a preferred embodiment of a method of processing an image according 

10 to the invention, filtering of the current detail image includes computation of weighted 

averages of pixel- values of the current detail image. Filtering of pixel-values in the current 
detail image includes computation of local weighted averages of pixel-values of the current 
detail image. In particular, filtering of the current detail image includes computation of a 
weighted average of neighbouring pixels in the current detail image with weights obtained by 

IS finding the direction in which the neighbouring pixel is located in the current detail image 
with respect to the pixel to be filtered, calculating the magnitude of the difference between 
pixel values in said direction at at least one of the coarser resolution levels. The higher this 
difference is, the lower the selected weight is. The weights have a low magnitude for pixel- 
values which are far apart in the current detail image, so that local averaging in the current 

20 detail image is achieved. As differences between, especially gradients of, pixel-values in the 
subsequent detail image and corresponding pixel-values in the current detail image are larger, 
the weights are smaller, so that in the current detail image filtering more predominantly 
retains the pixel-values of the current detail image. The weights are such that as differences 
between pixel-values in the subsequent detail image are larger, the computation of weighted 

25 averages in the current detail image takes pixel-values in a smaller region into account. 

Preferably, the weights have a smaller magnitude for relevant positions in the current detail 
image that are far apart, and the weights have a smaller magnitude as the gradient in the 
subsequent detail image has a larger magnitude at positions corresponding to said relevant 
positions in the current detail image. The weights also depend, in particular, on the relative 

30 direction of the gradient with respect to the difference between relevant positions in the 
subsequent detail image. Accordingly, the weights take directional structure in the image 
information in the subsequent image into account. Notably, the weights have lower 
magnimde as the gradient of the pixel-value is transverse to the difference between relevant 
positions in the subsequent detail image. To reduce noise in the image, in separate detail 
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images pixel-values are replaced by a weighted sum of the pixel value at issue in the detail 
image and the corresponding values for neighbouring pixels. If the vector pointing from the 
pixel at issue towards one of its neighbouring pixels is substantially parallel to the 
corresponding gradient at one of the lower resolution levels, the weight is low in order not to 

5 blur the image; if this vector points perpendicularly towards the gradient vector the weight is 
high to reduce noise. Thus, it is achieved to remove noise from the current detail image ani 
retain relevant image information, especially directional structure, as pixel-values are more 
locally smoothed, viz. averaged, in the current detail image, because, in this situation there 
is less corresponding structure in the subsequent detail image. Because directional structure 

10 in the image is appropriately taken into account, the method of processing an image 

according to the invention is particularly suitable for processing medical diagnostic images 
which image a patient's vascular system or which image a thin catheter which is inserted into 
a blood vessel. 

The weights take into account that the gradient of pixel-values in the 

15 subsequent detail image can vary from pixel to pixel. Thus, especially curved edges in the 
image are appropriately taken into account for controlling filtering. Edges in the image that 
represent important information may be curved; in medical diagnostic images that is often the 
case. Preferably, the magnitude of the difference between pixel values at at least one of the 
coarser resolution levels is calculated both for the pixel to be filtered and for the 

20 neighbouring pixel used for filtering. The weight is preferably chosen to be low if the 
magnitude of at least one of the two differences is high. For such curved structures, the 
direction of the gradient changes along that structure. Therefore, the gradient at the pixel to 
be filtered and at the neighbouring pixel may point into different directions. This means that 
a neighbouring pixel, being located perpendicularly to the gradient at the pixel to be filtered, 

25 may be located essentially in the direction of the gradient at that neighbouring pixel. 

Secondly, if the directional structure is not an edge but a line, the gradient 
magnitude may reach a minimum exactly on that line and may point along that line. Pixels 
near that line show a strong gradient perpendicular to that line. 

In each of the two cases, the weight for computation of a local average is 

30 reduced when the gradient - calculated from one of the lower resolution levels - at the pixel 
to be averaged or the neighbouring pixel is substantially in the direction of the vector 
pointing from the pixel to be filtered toward the neighbouring pixel. 

In a preferred embodiment of a method of processing an image according 
to the invention, the weights depend mainly on the angle between said gradient and a 
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difference between said positions 

When the gradient in the subsequent detail image has a large magnitude 
with respect to the noise level of the image, the weights are not very substantially dependent 
on the direction of the gradient relative to the difference between relevant positions in the 
subsequent detail image. In particular, it is avoided that almost all weights are close to zero 
when gradients have a direction which is not precisely perpendicular to the difference 
between relevant positions in the subsequent detail image. In this way. spurious sensitivity on 
the matrix structure of the subsequent detail image is avoided. 

In a preferred embodiment of a method of processing an image according 
to the invention, the coarsest resolution level is excluded from controlling the filtering. The 
gradient of the detail image at the lowest, i.e. coarsest, resolution level represents mainly 
linear ramps in the image, i.e. large areas with linearly increasing pixel-values in 
substantially one direction in the image. Thus, as the coarsest resolution level is excluded 
from controlling the filtering, influencing the filtering by linear ramps in the image is 
avoided. The sum of gradients in particular represents the position and strength of edges in 
the image. 

In a preferred embodiment of a method of processing an image according 
to the invention, the noise level of the image is derived mainly from image information in 
the image. When the noise level varies with the pixel-values in the image, during processing, 
notably filtering, automatically the actual noise level is taken into account. Deriving the noise 
level from the image information is relatively simple and requires uncomplicated 
computations. Often, there is a predetermined relation between the pixel-values and the noise 
level, e.g. due to a particular stochastic process which causes the noise. In particular for x- 
ray images, the noise level depends on pixel-values in the image. This property of the noise 
level is due to the quantum shot noise in the x-radiation employed to form the x-ray image. 
From the international application WO 96/01426 it is known per se that the noise level of the 
x-ray image can be derived, in essence, only from image information in the x-ray image. 
Often, a more accurate result for the noise level is obtained when, in addition, a few 
parameters pertaining to the setting of an x-ray examination apparatus employed to generate 
the x-ray image are taken into account. 

An image processor according to the invention comprises 

a decomposition-unit for decomposing the image into detail images at successive 
resolution levels, 

a filter for filtering detail images to generate filtered detail images. 
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- the filter being adjustable in dependence on at least one subsequent 
detailimage, said subsequent detail image having a resolution level which is lower than the 
resolution level of the current detail image, and 

a synthesizer for reconstructing a processed image from the filtered detail 

5 images. 

An image processor according to the invention is particularly suitable for 
performing the method of processing an image according to the invention. 

It is noted that the method of processing an image is particularly 
advantageously employed to process images that are generated by x-ray imaging, including 
10 computed-tomography, but also to process images generated by magnetic resonance imaging, 
ultrasound or electron microscopy. An x-ray examination apparatus comprises an x-ray 
detector, notably an image intensifier-television chain or an image sensor matrix. The x-ray 
detector derives an image signal, notably an electronic video signal, from the x-ray image. 
The image processor processes the image signal so as to generate a processed image signal. 
15 The processed image signal is suitable to display the image information in the x-ray image 
with a high diagnostic quality. 

It is noted that the functions of the image processor according to the 
invention may be carried out by a suitably programmed computer. Alternatively, the image 
processor may be provided with a special-purpose microprocessor that is designed to carry 
20 out the functions of the image processor. 

These and other aspects of the invention will be elucidated with reference 
to the embodiments described hereinafter and with reference to the accompanying drawing 
wherein 

Figure 1 shows a diagrammatic representation of an embodiment of an 
25 image processor according to the invention. 

Figure 2 shows a diagrammatic representation of another embodiment of 
an image processor according to the invention, 

Figure 3 shows a diagrammatic representation of a decomposition unit of 
the image processor according to the invention, 
30 Figure 4 shows a diagrammatic representation of a synthesizer of the 

image processor according to the invention, and 

Figure 5 shows a schematic representation of an x-ray examination 
apparatus wherein the invention is employed. 
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Figure 1 is a schematic diagram of an image processor according to the 
invention. An input image signal Qq is supplied to the decomposition unit 2 which generates a 
plurality of detail-image signals, notably frequency-band image signals at several resolution 
levels. The respective detail images contain image details at different spatial scales, and at 

5 lower resolution levels, and they contain coarser details. The signal levels of the input image 
signal Qq represent pixel-values of the input linage. The operation of the decomposition unit 
is further elaborated with reference to Figure 3. In the example, four decomposition levels 
are shown, but in practice more resolution levels may be employed, for example for an input 
image having Ik^ (1024 x 1024) pixels, five or six resolution levels may be used. At each 

10 resolution level, detail-image signals, notably a high-pass frequency band signal Ajt, are 

formed and at the lowest resolution level, also a low-pass frequency -band signal is formed. 
Increasing index k indicates lower resolution levels. At the highest resolution level (i.e. the 
0th level) the high-pass frequency-band image signal Hq is applied to a filter 41 which 
supplies the filtered detail image, notably a filtered high-pass frequency band signal, to the 

15 synthesizer unit 3, At lower resolution levels, the high-pass frequency-band image signals hj, 
are applied to filters 42,43 which supply the filtered detail images to the synthesizer unit 3. 
Also, the low-pass frequency-band image signal /j of the lowest resolution level is supplied to 
the synthesizer unit 3. The processed image signal Qq is reconstructed by the synthesizer unit 
3 from the various filtered detail image signals at the respective resolution levels. The 

20 operation of the synthesizer unit 3 is further elaborated with reference to Figure 4. 

The filters 41, 42, 43 compute the local weighted average at respective 
resolution levels of pixel-values of the respective detail-image signals h0,h„h2 and hj 
respectively. The respective filters 41,42,and 43 which generate filtered detail images hf, are 
adjusted on the basis of gradients gj^+j in the detail image of the next coarser resolution level. 

25 

Notably the pixel-values of the filtered detail-image h^' are computed as 
the weighted average 

hf. ix) =J2a {Ax, X, at^i) (x+Ax) 

Here 5 is a usually small region around the position jc, preferably S consists of the 
neighbouring positions of the current position x in the current detail image. The weight 
30 factors a are normalised so that their sum over the region S equals unity. The weight factors 
a are decreasing functions of the local gradient and the distance vector Ax, a particular form 
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of the weight factors is 



a (Ax, X, sr^^i) =/S (Ax) r (x) , Ax] r [g^^^ (x+Ax) , Ax] 



where the smoothly decreasing function r is, for example, a Gaussian function 



r(5r. Ax) =exp 



r, (gr-Ax)^ l 
L ^(llsrll) J 



The rate of decrease of the function r is determined by the denominator function v. The 
denominator function increases with increasing gradients, so that the averaging effect is 

5 stronger more locally in the current detail image when there is a corresponding large gradient 
in the next coarser resolution level. Moreover, it is advantageous to take into account that 
there are usually contributions due to noise in the gradients. Notably, in an x-ray image, 
quantum shot noise gives rise to a noise component in the input image signals that is 
conveyed to the detail-image signals. The denominator function is arranged such that for very 

10 large gradients averaging is effective over a wide region beyond the position at issue, and 
when the gradient is below the noise level, averaging is effectively concentrated on a close 
neighbourhood of the position at issue. It appears that a suitable denominator function is 



Here a(g)^ is the variance of the gradient which serves as a fair estimate of the noise 
component in the gradient. The scalar parameters and L are adjustable so as the set the 
15 relative magnitude of the various terms in the denominator function. The scalar parameters 
can be set empirically so that filtering has a desired directional sensitivity in order to retain 
elongate structures in the image The directional sensitivity of the denominator function is 
controlled by the parameter L, i.e. for very large gradients the denominator function depends 
mainly asymptotically on the angle between the gradient and the vector Ax, i.e. 



Figure 3 is a schematic diagram of the decomposition unit 2 of the image 
processor of Figure 1. The image signal Qq is supplied to low-pass down-sampling filters 101 
and 102 which perform low-pass filtermg in two directions in the image (denoted as x and y) 



M||gr||)= c-t 
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and a down-sampling operation, for example by a factor of 2. Down-sampling is simply 
performed by omitting every second pixel-value. The low-pass down-sampled signal is 
supplied to an interpolator 200, The interpolator 200 is formed here by a pair of low-pass 
up-sampling filters 111,112 which also perform an up-sampling operation in both x and y 

5 directions. The operation as an interpolator of the low-pass up-sampling filters 111,112 
serves to insert zeros between successive pixel-values in the signal Uj and to perform a 
smoothing operation. The low-pass up-sampling filters 101,102 and low-pass down-sampling 
filters 111,112 may have the same or different cut-off frequencies. The interpolator 200 
derives a low-pass frequency band image signal Z^, which is subtracted by subtracter 120 

10 from the input image signal Qq so as to form the high-pass frequency-band image signal /iq. 
The signals and Iq include image information having variations on at least a spatial scale 
that corresponds to the cut-off frequency of the low-pass down-sampling filters 101 and 102. 

The signal is subsequently passed to a pair of low-pass down-sampling 
filters 103 and 104 to form a signal A low-pass frequency band image signal 1 is derived 

15 from the signal ^2 by means of an interpolator 201 which comprises a pair of low-pass up- 
sampling filters 113,114. The signals ^2 h include image information having variations 
on at least a spatial scale that corresponds to the cut-off frequency of the low-pass down- 
sampling filters 103 and 104. By means of subtracter 121 the high-pass frequency band 
image signal for the 1st resolution level is derived from the signals and /y, viz. 

20 hi = arh^ 

The signal ^2 is subsequently passed to a pair of low-pass down-sampling 
filters 105 and 106 to form a signal aj. A low-pass frequency band image signal is derived 
from the signal Oj by. means of an interpolator 202 which comprises a pair of low-pass up- 
sampling filters 115,116. By means of subtracter 122 the high-pass frequency band image 

25 signal A2 for the 2nd resolution level is derived from the signals and /j. The signals Cj and 
I2 include image information having variations on at least a spatial scale that corresponds to 
the cut-off frequency of the low-pass down-sampling filters 105 and 106. 

The signal is subsequently passed to a pair of low-pass down-sampling 
filters 107 and 108 to form a signal a^, A low-pass frequency band image signal Ij is derived 

30 from the signal by means of an interpolator 203 which comprises a pair of low-pass up- 
sampling fillers 117,118. By means of subtracter 123 the high-pass frequency band image 
signal hj for the 3rd resolution level is derived from the signals and Zj. The signals and 
Ij include image information having variations on at least a spatial scale that corresponds to 




wo 98/55916 



PCT/IB98/00850 



11 

the cut-off frequency of the low-pass down-sampling filters 107 and 108. 

It will be apparent to the skilled person that the decomposition into 
frequency-band image signals for successive resolution levels can be continued beyond the 
four levels (Ost to 3rd levels) shown in the embodiment of Figure 3. 
5 Figure 4 is a schematic diagram of the synthesizer 5 of the image 

processor of Figure 1. The synthesizer 5 comprises an adder 129 arranged to add the low- 
pass frequency-band image signal Ij of the lowest (in this case the 2nd) resolution level to the 
processed high-pass frequency-band image signal A j so as to form a signal af. An 
interpolator 149 formed by low-pass up-sampling filters 139,140, derives a low-pass 
10 processed frequency-band image signal of the 1st resolution level. The synthesizer 5 

further comprises an adder 130 arranged to add the low-pass frequency -band image signal /^^ 
of the lowest but one (in this case the 2nd) resolution level to the processed high-pass 
frequency-band image signal /j2~ so as to form a signal ^2". An interpolator 150 formed by 
low-pass up-sampling filters 141,142, derives a low-pass processed frequency-band image 
15 signal /y^ of the 1st resolution level. Another adder 131 adds the processed high-pass 

frequency-band image signal of the 1st resolution level hj" to the signal //" so as to form a 
signal Of". A low-pass processed frequency-band image signal of the 0th resolution level lo" 
is derived by a further interpolator 151 from the signal Qj". Finally, the filtered output 
image signal Gq'^ of the 0th resolution level is formed by still another adder 132 which adds 
20 up the signals Iq" and h'o- The signal levels, viz. the signal amplitudes of the filtered output 
image signal a^" , represent pixel- values of the filtered output image. The interpolator 151 
comprises two low-pass up-sampling filters 143,144. The interpolators 150 and 151 operate 
in the same manner as the interpolators 200. 201 of the decomposition unit 2. It is noted that 
the frequency characteristics of the low-pass up-sampling filters of the interpolators 150,151 
25 should be accurately matched at separate resolution levels to the frequency characteristics of 
low-pass up-sampling filters 111 to 118 of the respective corresponding resolution levels. 
Said accurate matching is required so as to avoid loss of image information in the combined 
decomposition-synthesisation process. 



30 provided with an image processor according to the invention. The x-ray examination 

apparatus 11 comprises an x-ray source 12 for irradiating an object 13, for example a patient 
who is to be radiologically examined, with an x-ray beam 14. Because of local variations of 
the x-ray absorption in the patient, an x-ray shadow image is formed on an x-ray sensitive 



Figure 5 is a schematic diagram of an x-ray examination apparatus 
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face 15 of an x-ray detector 16. An electronic image signal is derived from the x-ray image 
by the x-ray detector. In particular the x-ray detector comprises an x-ray image intensifier 17 
having an entrance section 18 which incorporates the x-ray sensitive face 15. The x-radiation 
incident on the entrance section 18 is converted into an image carrying electron beam which 
generates a light-optical image on the exit window 19, The entrance section includes the x- 
ray sensitive face, having the form of a scintillator screen 30, and a photocathode 31. The 
incident x-rays are converted in the scintillator screen into e.g. blue or ultraviolet light for 
which the photocathode 31 is sensitive. The electron beam 32 generated by the photocathode 
31 is imaged by an electron-optical system 32 onto a phosphor screen 33 that is disposed on 
the exit window 19. The electron optical system includes the photocathode 31, a hollow 
anode 34 and a number of electrodes 35! A camera 20 is arranged to pick-up the light-optical 
image; to that end an optical coupling 21, e.g. a lens system, is provided to optically couple 
the camera 21 to the exit window 19. The electronic image signal at the output of the 
camera is supplied as the input image signal to the image processor 1 according to the 
invention. The output signal a^o represents a filtered output image in which noise is 
substantially reduced while small details are preserved. In practice, a noise reduction of 
4-5dB can be achieved. Therefore, the filtered output image has a high diagnostic quality 
even when a low x-ray dose is employed. 

Reference is made again to Figure 1 for a discussion of the adjustment of 
the filters at various resolution levels, Subtractors 51,52,53 are provided which derive local 
differences d^.^, from filtered detail image signals at the next coarser resolution levels. For 
example, the subtracter 51 together with an addition unit 54, an up-sampler 56 and a low- 
pass filter 59 derives a local gradient which is employed to control the filter 41. To that 
end, the low-pass filter 59 is coupled to a control input 61 of the filter 41. More particularly, 
the subtracter computes a difference between pixel- values, the difference < generated by the 
subtracter 51 is interpolated with the corresponding gradient g^^t from the next coarser 
resolution level, by way of addition by the addition unit 54 and by employing the up-sampler 
56 which inserts zeros in the signal representing the gradient gf,. The up-sampling is carried- 
out to render the signal representing the gradient suitable to be combined with the filtered 
detail image signal at the next lower (less coarse) resolution level. Further, a low-pass filter 
59 smooths the up-sampled signal so as to achieve a gradual interpolation of the gradients of 
adjoining resolution levels and divides the gradient by a factor so as to compensate for a 
difference between effective distances between pixels in different detail images at respective 
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resolution levels. Notably, the addition unit 54, the up-sampler 56 and the low-pass filter 59 
form an interpolator. Similarly the filter 42 of the 1st resolution level is adjusted on the basis 
of the gradient at the 2nd resolution level which is derived by a subtracter 52 together 
with an addition unit 55, an up-sampler 57 and a low-pass filter 60. Further, the filter 43 at 

5 the 2nd resolution level, i.e. the highest but one (coarsest but one) resolution level is 
adjusted on the basis of the gradient gj at the highest resolution level which is not 
interpolated. As an alternative, the local gradients gj 2 may be derived from the high-pass 
frequency-band image signals /z; 2 as indicated by the dashed lines in Figure 1. 

Figure 2 shows an alternative embodiment of a signal processor according 

10 to the invention. The subtracters 51,52,53 and the interpolation units are arranged so as to 
compute gradients of pixel-values of interpolations of filtered detail images at successive 
resolution levels. Thus, in fact the interpolation units and the subtracters are interchanged as 
compared to the embodiment of Figure 1. At the lowest level, the gradient of the detail 
image hj is computed by a separate subtracter 65. The gradient of the detail image at the 

15 coarsest resolution level is used to control the filter 44 for the coarsest resolution level. 



wo .98/55916 



14 



PCT/IB98/00850 



Claims 



1 . A method of processing an image comprising the ^steps of 
decomposing the image into detail images at successive resolution levels, 
filtering detail images to generate filtered detail images, 

- filtering a current detail image in dependence on at least one subsequent detail 
5 image, said subsequent detail image having a resolution level which is lower than the 

resolution level of the current detail image, and 

reconstructing a processed image from the filtered detail images. 

2. A method of processing an image as claimed in Claim 1 , wherein 
filtering of the current detail image takes place in dependence on a next detail 

10 image having a highest resolution level which is lower than the resolution level of the current 
detail image. 

3. A method of processing an image as claimed in Claim 1, 

wherein filtering of the current detail image depends on differences between 
pixel-values of at least one subsequent detail image having a resolution level which is lower 
15 than the resolution level of the current detail image. 

4. A method of processing an image as claimed in Claim 3, wherein 
filtering of the current detail image includes computation of weighted averages 

of pixel-values of the current detail image, with weights depending on 

- differences between pixel-values of pixels of said subsequent detail image, and 
20 - differences between positions of said pixels in said subsequent detail image. 

5. A method of processing an image as claimed in Claim 4, wherein 

the weights depend on gradients of pixel-values of said subsequent detail image, 
said gradients relating to several positions in said subsequent detail image. 

6. A method of processing an image as claimed in Claim 4, wherein 
25 - the weights are computed as a decreasing function of 

- differences between pixel-values of pixels of said subsequent detail image and 

of 

- differences between positions of said pixels in said subsequent detail image. 

7. A method of processing an image as claimed in Claim 4, wherein 
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for a gradient of pixel-values at respective positions in said subsequent detail 
image, said gradient being substantially larger than a threshold that relates to a noise level of 
the image, the weights depend mainly on the angle between said gradient and a difference 
between said positions 

5 

8. A method of processing an image as claimed in Claim 1 or 5, wherein 

filtering of the current detail image takes place in dependence on a sum of 

gradients of pixel-values of detail images from a group of detail images not including the 

detail image at the coarsest resolution level. 
10 9. A method of processing an image as claimed in Claim 1 wherein 

the noise level of the image is derived mainly from image information in the 

image. 

10. An image processor comprising 

a decomposition-unit (2) for decomposing the image into detail images at 
15 successive resolution levels, 

a filter (41, 42, 43) for filtering detail images to generate filtered detail images, 
- the filter being adjustable in dependence on at least one subsequent detail 
image, said subsequent detail image having a resolution level which is lower than the 
resolution level of the current detail image, and 
20 - a synthesizer (3) for reconstructing a processed image from the filtered detail 

images. 

11. An x-ray examination apparatus comprising an image processor as claimed 
in Claim 10. 
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